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Abstract 



The mean sky-averaged Comptonization parameter, y , describing the scattering of the CMB by 
hot gas in clusters of galaxies is calculated in an array of flat and open cosmological and dark 
matter models. The models are globally normalized to fit cluster X-ray data, and intracluster 
gas is assumed to have evolved in a manner consistent with current observations. We predict 
values of y lower than the COBE/FIRAS upper limit. The corresponding values of the overall 
optical thickness to Compton scattering are ^ 10~ 4 for relevant parameter values. Of more 
practical importance are number counts of clusters across which a net flux (with respect to the 
CMB) higher than some limiting value can be detected. Such number counts are specifically 
predicted for the COBRAS /SAMBA and BOOMERANG missions. 

1 Introduction 

Compton scattering of the cosmic microwave background (CMB) radiation by hot gas in clusters 
of galaxies - the Sunyaev-Zeldovich (1970; hereafter S-Z) effect - affects both the spectrum 
and the angular distribution of the radiation. The spectral change of the CMB across each 
cluster translates into a superposed average change of the spectrum across the sky. When the 
intensity (or temperature) is differentially measured, the radiation will appear anisotropic, on 
characteristic angular scales of several arcminutes, reflecting the typical sizes of clusters that 
contribute most to the anisotropy. This cluster-induced anisotropy is an important component 
of the full anisotropy on such scales, as has been shown in many studies (see, e.g. , Rephaeli 
1995a, and references therein). 

Less studied is the overall spectral effect of clusters on the CMB. Whereas the anisotropy is a 
measure of the rms value of the Comptonization parameter y [defined in Equation (12) below], 
the overall spectral effect results from the superposed effects of individual clusters over the sky. 
This parameter is of prime interest in the characterization of spectral deviations from a pure 
blackbody spectrum. Currently, the COBE/FIRAS database provides the best measure of such 
non-Planckian distortions by constraining the degree of Comptonization to be y < 1.5 x 1CT 5 at 
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the 90% statistical significance level (Fixsen et al. 1996). The significance of this value stems 
from the fact that it limits the cumulative effects of energy release processes in the early universe 
(at z ^ 3 x 10 6 ), and the superposed effects of hot intergalactic (IG) and intracluster (IC) gas 
(Wright et al. 1994). In order to predict the contribution to the Comptonization parameter due 
to a population of evolving clusters, some assumptions need to be made about the evolution of 
groups and clusters (which, for brevity, will hitherto be collectively referred to as clusters) of 
galaxies, and basic characteristics of their hot gas. Under these assumptions the distribution 
of the Comptonization parameter in clusters at different redshifts and in a given mass range 
can be calculated. From this, the mean Comptonization due to clusters can be reasonably well 
estimated. 

Here we present the results of an investigation of the effects on the CMB of a population of 
evolving clusters, described in the background of various flat and open cosmological models. In 
a previous paper (Colafrancesco et al. 1994) we have reported results of a calculation of the 
predicted rms fluctuations in the CMB induced by gas in clusters in a set of flat cosmological 
models. Because our basic approach here is similar to that detailed in the latter paper, our 
discussion - in §2 - of the general background and modeling of clusters and of their IC gas will 
be brief. The basic definitions of the sky-averaged quantities we calculate here are given in §3. 
We also calculate the expected cluster number counts in the mm and sub-mm bands. 

Long duration balloon flights and space missions dedicated to the study of CMB anisotropy 
are expected to reach sensitivity levels of ~ 10 /iK per pixel. The COBRAS/SAMBA mission 
(hereafter C/S; see e.g. Mandolesi et al. 1995), has been selected as the next Medium-Sized 
Mission of the ESA Horizon 2000 Program. The eight channels of this experiment cover a 
wide frequency range, from 35 to 714 GHz. In particular, the 140, 222 and 400 GHz channels 
of the High Frequency Instrument (HFI) are particularly relevant for measurement of the S-Z 
effect in clusters. The first frequency is on the Rayleigh- Jeans (R-J) side, where the intensity 
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change (across a cluster) is negative; the second is near the crossover frequency (where the 
thermal effect vanishes), and the third is on the Wien side, where the effect is positive (for 
a discussion of the significance of high-frequency measurements, see Rephaeli, 1995a). The 
corresponding angular resolutions at these frequencies are 10.5', 7.5', 4.5'(FWHM), respectively. 
A long duration balloon-borne experiment BOOMERANG (see Lange et al. 1995) is expected 
to cover a 10° x 10° patch of the sky with a FWHM resolution of 5'at 150 GHz. The projected 
sensitivity levels and angular resolution of these experiments are very suitable for a survey of 
a large number of clusters. In §4 we present the results of our calculations giving the predicted 
ranges of values of basic quantities which can be determined in light of these observational 
capabilities. Finally, in §5 we assess these results and summarize our main conclusions. 



2 Cosmological Evolution of Clusters 

2.1 Mass and Redshift Distribution of Clusters 

We adopt the simple spherical collapse picture for cluster formation, according to which a 
homogeneous, spherical perturbation detaches from the Hubble flow at time t m , collapses at 
time t c ~ 2t m , and virializes at time t v ~ 3t m . The relative density contrast at an initial 
redshift z iy Si jV , depends on the cosmological model. Under the assumption of linear growth, 
the density contrast at t v is 5 V = 8 ijV D(t v )/D(ti), where D(t v ) is the linear growth factor in the 
chosen cosmology. The relevant formulae for low density (open, or vacuum dominated) models 
are given in the Appendix. For Q 1, 5 V tends to the standard value of 2.2, independent of 
t v . The actual, non-linear density contrast is A = p/p, where p is the perturbation density and 
p is the background density at the time of virialization. In the Appendix we also briefly review 
how to derive A in low density cosmological models. For z v — > oo, A, evaluate at z v tends to 
the standard value ~ 400, found in a flat universe. 
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The mass and redshift distribution of clusters, N(M,z), can be determined from the Press 
& Schechter (1974; hereafter P&S) theory: 



where p is the background density at redshift z, M is the total cluster mass, a(M, z) is the 
rms of the linear density fluctuation field at z, smoothed over the region containing M, and 
S v is the linear density contrast of a perturbation that virializes at z. The variance of density 
fluctuations of mass M is given in the standard relation: 



a 2 {R,z) = D 2 {z) J dink k 3 P{k) 



(2) 



kR 

with M = 47rpi? 3 /3. We use the fitting formulae for the power spectra, P(k), given in Holtzman 
(1989) in order to calculate the mass variance in different cosmological models. Here we consider 
flat cold (CDM) and mixed (MDM) dark matter models, low density CDM models with a 
cosmological constant (CDM+A), as well as open, pure baryonic models with isocurvature 
initial conditions (BDM). 

If the power spectrum of density fluctuations is a power-law of the spatial frequency [i.e. 
P{k) = Ak n ], then a(M, z) = (l/b)(M/M Q )~ a D(z), where M is the mass contained in a sphere 
of 8/i _1 Mpc radius, b is the biasing factor, and a = (n + 3)/6. For this power-law spectrum, 
the mass distribution assumes the well known form: 



J p n + 3 S v b fM\ 



a-2 



1 b 2 v b 2 (M\ 2a 



2D 2 (z)\M J J' (3) 

where 1 < J < 2 takes into account possible secondary infall of mass into the cluster after its 
initial collapse and virialization. 

The fact that N(M, z) depends on the product S v b, and not separately on 5 V and b, somewhat 
simplifies the fit to X-ray data. We have determined the values of X and 5 v b in various dark 

matter models by fitting to the observed cluster X-ray luminosity function of Kowalski et al. 
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(1984). The results of our analysis (see Colafrancesco & Vittorio 1994 for details) are showns in 
Table 1. Using a P&S theory where X and 5 v b are fitted to the X-ray luminosity (or temperature) 
functions, weakens the predictive power of the theoretical models based on the linear theory. 
However, the fitted values of X and 5 v b should contain the relevant informations of a more 
realistic cluster formation picture. 

2.2 Properties and Evolution of Intracluster Gas 

There are various open issues pertaining to the formation of the hot gaseous cores of clusters. 
For our purposes here we can adopt the following simplified approach: shortly after a cluster 
forms and virializes, a gaseous core forms (probably as a result of tidal galactic interactions and 
other gas stripping processes) with the hot gas in hydrostatic equilibrium in the potential well 
of the cluster. To avoid introducing a large number of free parameters, we will simply scale the 
gas properties to those of the cluster. This is suggested by X-ray observations of local clusters 
(see, e.g. , Jones & Forman 1992). The gas mass is taken to be a fraction / of the total cluster 
mass, and the gas density profile is assumed to have the commonly adopted form 

n(r)=n [l + (r/r c y]-W\ (4) 

where n is the central electron density and r c is a core radius. The observed values of (5 range 
roughly from 0.5 to 0.7. Here we use the value (3 — 2/3 which is particularly convenient in 
analytic calculations. 

The radial extent of a cluster is taken as R = pr c . The mass within the outer radius R(p) is 

M(p) = 3M [p-tg- 1 p}, (5) 
where M = (47r/3)r^p 0) an d Po is the central total mass density of the cluster. 

Assuming that cluster collapse is self-similar, we can infer the mass and redshift dependence 
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of r c from the scaling law: r c 



3M/4np b A 



r c (Q ,M,z) = 



1.29 h' 1 Mpc 
P 



1/3 

(1 + z) . This gives 
M A(ft o = l,z = 0) 



1/3 ]_ 



1 + 2 



(6) 



Hereafter we fix p = 10 to recover values of the IC gas core radii consistent with observations 
(see, e.g. , Sarazin 1988, Jones & Forman 1991). In particular, for a local cluster of 10 15 h^ 1 M Q 
we get r c = 0.12 and 0.16 hr 1 Mpc for Q — 1 and f2 = 0.2, respectively. 



The gas is assumed isothermal at the virial temperature T oc Mj R. This gives 



T = 8.7-10 7 (l + z)( 



M 



V$*>hr x M, 



v2/3 







^qA(^o^) 

A(n = i,z = o) 



1/3 



K. 



(7) 



Although it is known that the gas mass fraction depends on z and M, little is currently known 
on the exact form of these dependences. We adopt the simple parametrization (described in 
detail in Colafrancesco & Vittorio 1994; see also Cavaliere, Colafrancesco and Menci 1993) 
which is based on the results of analyses of the Einstein Medium Sensitivity Survey (EMSS) 
data (Gioia et al. 1990; Henry et al. 1992) which seem to indicate a decrease in the number 
of bright clusters with redshift, and an analysis of a local cluster sample (David et al. 1990): 



f = fo 



M 



V t \Z 



(8) 



.io 1 ^- 1 m ; \t 

Here t Q is the age of the universe, and the normalization to f Q ~ 0.1, is based on a local, rich 
cluster sample. Values of rj and £ are listed in Table 1. 



3 Cluster Comptonization 

In the non-relativistic limit, the effect of scattering of the CMB by hot gas depends linearly on 
the cluster Comptonization parameter: 

Vc = / nM ' ( 9 ) 

7Tl e C J irnin 
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where ot is the Thomson cross section, and the integral is over a line of sight through the cluster. 
In the exact relativistic treatment, the dependence of the effect is not exactly linear in y. For 
low values of the optical thickness r to Compton scattering, the dependence on r is linear, but 
the dependence on T is more complicated, even for the observed range of IC gas temperatures, 
roughly 3 -r- 15 keV. The higher the gas temperature and the observing frequency, the larger 
is the deviation of the intensity change from the non- relativistic value (Rephaeli 1995b). This 
must be taken into account in the analysis of high frequency observations of individual rich 
clusters (Rephaeli 1995a, Holzapfel et al. 1996). Here we are interested in the integrated 
effect due to many clusters, and since we expect this to be largely dominated by the numerous 
low richness (and, correspondingly, low temperature) clusters (Colafrancesco et al. 1994), we 
retain here the considerable degree of simplicity which results from expressing the full spatial 
dependence linearly in the Comptonization parameter. The overall effect of this simplification 
on our results is assessed in the Discussion. 

The cumulative Comptonization parameter 2/(7) along the line of sight (los) 7 is the sum 
over all clusters whose gaseous spheres are intersected by this los. Identifying a cluster in the 
ensemble by its mass M m and redshift zi, the expression for 2/(7) can be written as (Cole and 
Kaiser 1988): 

V(l) =^Zn hm y (M m ,z l )Ci\i -ii\,M m ,zi) . (10) 

l,m 

Here n^ m is the occupation number of clusters in the M — z space, while 

y {M m ,zi) = (2tg~ 1 p)(kT/m e c 2 )(T T n r c - (11) 

is the value of the Comptonization parameter along a los, 7/, through the center of the cluster 
(with mass M m at redshift z{). The angular profile (, obtained by integrating the truncated 
profile of Equation (6) along different directions is 



<Z(\7-7 l \,M,z)=tg- 1 



8 



where cos(0) = 7 • % 9 < p9 c , and 9 C = 9 C (M, z) is the angle subtended by r c . 

The distribution of IC gas has characteristic angular scales reflecting the apparent sizes 
of clusters at various redshifts (unlike the case of a possible uniform intergalactic medium). 
Nevertheless, even though the distribution of values of the Comptonization parameter depends 
on the beam size, its true sky-averaged value is beam- independent. To find the mean value of 
the Comptonization parameter, we average 2/(7) first over the sky and then over the ensemble 
of cluster distributions. The sky and the ensemble average act on the profile and the phase 
space density respectively, yielding 

ZU^z) = ^r p {p-tg- 1 p). (13) 

Altogether, the mean Comptonization parameter averaged over the ensemble of clusters and 
over the sky is 

r dV r 

V = J ~^dz J N(M, z)y (M, z)((M, z)dM . (14) 

In Table 2 we list values of y, in the selected open and flat models with the parameters of 
Table 1 and M min = 10 13 M . The values of y are generally well below the FIRAS upper limit, 
1.5 x 1(T 5 (Fixsen et al. 1996). 

Because values of most of the parameters are poorly known, it is important to consider 
the full parameter range. In Table 2 we also list predicted values of y obtained with no gas 
evolution {i.e. rj = 0, £ = 0), and with different choices for M min . Based on these estimates 
we conclude that the assumed degree of IC gas evolution reduces y by a factor of few, except 
in the BDM models where this factor can be very significant. The decrease in the degree of 
Comptonization with increasing M min is dramatic only in the BDM Q = 0.1 model. The 
predicted value of y in the BDM Q = 0.1 model with no gas evolution is inconsistent with the 
current FIRAS limit, even if M min is increased by a factor of ~ 3 to M rnin = 1O 13 M . 
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As shown in Figure 1, most of the contribution to the mean Comptonization comes from 
clusters at z ^ 2, if there is gas evolution. For open CDM models without gas evolution there 
is a substantial contribution to y up to z ^ 5. 

Another useful quantity which can be readily calculated using the above formalism is the 
mean optical thickness to Compton scattering due to gas in clusters. Since the optical thickness 
of a cluster is r = y ■ m e c 2 /kT, the above expressions can be used to calculate T by simply 
replacing y (M,z) in Equation (15) with t q (M,z) = y (M,z)m e c 2 /kT, the optical thickness 
along a los through the center of a cluster of mass M at redshift z. Values of the mean sky- 
averaged optical thickness of gas in clusters are listed in Table 2. These give a measure of the 
minimal degree of scattering of the radiation from distant sources. The sensitivity to M min is 
larger in this case: nonetheless, values ^ 10~ 4 are obtained for most of the models. Only in 
the BDM models with n = — 1 and without gas evolution does r attain values of <: 10~ 3 . 

4 Cluster Counts 

With respect to the incident radiation field, the change of the CMB intensity across a cluster 
can be viewed as a net flux emanating from the cluster. The flux is negative below the crossover 
frequency and positive above this characteristic frequency (~ 217 GHz in the nonrelativistic 
limit). While the main emphasis so far has been the measurement of the S-Z effect in individual 
clusters, the capability to observe a large number of clusters in a satellite survey enhances 
interest in S-Z number counts. The number of clusters observable in sub-mm bands is of 
particular interest in view of forthcoming experiments dedicated to the study of the small scale 
structure of CMB anisotropy. In the nonrelativistic limit, the change of spectral intensity across 
a cluster of a given y, observed at a frequency x = hu/kT C MB [Tomb = (2.726 ± 0.010) K; 
Mather et al. 1994), is 

2(/cT C mb) 3 , v n _v 

= (he) 2 ' s( x )y , (15) 
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where 



4 T 

x e 



9( x ) = (e x _ 1)2 • bcotMx/2) - 4] • (16) 
The differential flux measured at a given frequency from the cluster is 

AF„(7,)=/ d^dT-Til.^A/^T), (17) 

J47T 

where i? s is the angular response of the receiver, say with a Gaussian beam of dispersion <r B , 
whose axis coincides with the los to the cluster center. The angular dependence of AI U is fully 
contained in the cluster profile, and so AF U can be computed by convolving this profile with 
the response of the receiver. From Equations (17) and (19) the beam convolved cluster signal 
is: 

AF„ = 2{1 ^™ 2 b)3 g(x)y (M, z)E(M, z) , (18) 

where 

Z = J dQR s (\^-^\)Cmi\,M,z). (19) 

The receiver measures the flux integrated over its passband E(v): 

_ J dvAF v E{y) 
V ~ fduE(u) ' (20) 

where E{y) is the frequency response of the instrument. Thus, the predicted number of clusters 



with a net flux > AF^ is 

N(> AF U ) = / — dz _ dMN(M, z) , (21) 
J dz Jm(af„,z) 

The lower bound of the mass integral, M, is determined from the requirement that the source 



flux is > AF V . Equation (20) provides the mass dependence of the flux, AF„(M), a monotonic 
function of the mass at a given z. 

In our specific estimates of the predicted number counts for the C/S experiment, we take 

the spectral response to be uniform over the two passbands centered on 140 and 400 GHz with 
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widths Av/v = Ax/x = 0.4 and 0.7, respectively. In Figures 3 and 4 we show the cluster 
number counts at the C/S's 140 and 400 GHz bands in the different models, with and without 
IC gas evolution. In Figure 5 we compare the counts expected in few models for the 400 GHz 
channel. 

The receivers of the C/S HFI and BOOMERANG are bolometers with a Noise Equivalent 
Power (NEP) of « 10~ 17 WHz~ 1/2 . The limiting flux is expected to be 

Avy/iAe' V ; 

where t is the integration time in seconds, A is the effective mirror area in squared meters, and 

e is the total (optical and electric) efficiency of the system. For C/S, NEP values of 1.0 x 10~ 16 

and 2.7 x 10~ 17 WHz~ l/2 are expected at 140 and 400 GHz, respectively. The C/S telescope has 

an effective diameter of 1 meter, and total efficiency of pa 0.30, so that with an integration time 

of one year (assuming full sky coverage) we have F v as 160 and 20 mJy at 140 and 400 GHz, 

respectively. A similar calculation of the BOOMERANG limiting flux - with the same values 

of A and e but with NEP « 2 • lO" 17 W Hz" 1 / 2 and t « 225 s per pixel - yields T v ° ise « 20 

mJy. In Tables 3 and 4 we give the number of clusters with flux greater than p™ lse anc l 3~p^, mse , 

corresponding to 1 and 3 sigma detection, respectively, as expected in the various models. Note 

that these numbers are obtained by requiring that the flux collected by the receiver from the 

center of the cluster is greater than the limiting flux. This estimate could be conservative for 

extended clusters if the flux emanating from the central region is under the detection limit. 

However, a further smoothing of the observed map could help in extracting additional clusters 

out of the noise (the noise of the smoothed map is lower by the factor smoot h/0 o b s , where 9 smoot h 

and 9 b s are the resolutions of the smoothed and observed maps, respectively). This is reflected 

in Figures 3 and 4, where we show the cluster number counts expected after smoothing the 

original map. The noise level of the smoothed maps at 400 GHz is 10.3, 3 and 1.5 mJy for final 

resolutions of 10, 30 and 60 arc-minutes, respectively. Correspondingly, the number of clusters 
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that can be detected increases. 



Even under the more conservative assumption of having just the central pixel above the 
limiting flux, the predicted counts are fairly large for C/S for which a full sky coverage is 
assumed; the quoted numbers for the BOOMERANG'S 150 GHz band (Ax/x = 0.2) refer to 
a 10° x 10° patch of the sky. Thus, given that the S-Z effect can be identified by virtue of its 
characteristic spectral signature, these experiments can, in principle, produce S-Z catalogs of 
clusters. Correlation analyses of S-Z and X-ray measurements of clusters will provide useful 
information on cluster evolution, and possibly also on cosmological scenarios. 

In Figure 6 we show the redshift distribution of the predicted number of clusters seen above 
the limiting flux at 400 and 140 GHz. For flat (CDM and MDM) models 99% of the predicted 
clusters have redshifts less than z ~ 0.2, quite independently of the degree of gas evolution. 
In low density (either open or with cosmological constant) models clusters are expected also 
at higher redshifts [z ^ 1). However, considering gas evolution in these low density models, 
~ 90% of the predicted clusters have z ^ 0.3. 

Most of the detected clusters are expected to be contained just in one pixel of the mm and 
sub-mm maps, and so it will be difficult to clearly identify clusters under these circumstances. It 
is therefore important to determine also the distribution of angular sizes of clusters whose fluxes 
are higher than the above limiting values. This is accomplished by calculating the distribution 

N(9 FWHM ,z) = N(M,z)^f-^— (23) 

where 9 C = Dat c is the angle subtended by the core radius of a cluster, which depends on M 
through Equation (8), Da is the angular diameter distance, and Ofwhm = 6fwhm{@c) is the 
FWHM of the beam convolved cluster profile: 

2 , n2- 



((,, ,„) = _L J mmexp { _ *_+».) ,„ (*) (24) 
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In Figure 7 we show the quantity 

POO r 

jV(> Qfwhm) = / d9' FWHM / dzN(9' FWHM , z) , 

JOfwhm J 

where the double integral is performed so that the mass of a cluster with a flux > p n ° lse j s 
larger than a minimum value M (see Equation 23); correspondingly, at a given z, such a cluster 
has a FWHM angular size larger than 9 FWHM - 

5 Discussion. 

The approach described in this paper is phenomenological, as it is based on the normalization 
of the predicted cluster number counts so they reproduce the locally observed distribution. We 
consider our normalization to the XRLF to be quite robust, in the sense that a statistically 
significant fit to the data is obtained over a range of X-ray luminosities from ~ 10 43 to ^ 10 45 
erg s _1 . The quality of this fit is low for MDM fi„ = 0.3 models with scale invariant initial 
conditions, unless n is increased to values ~ 1.2 -j- 1.4 (see Lucchin et al. 1996 for a discussion). 
BDM models with n in the range — l-j-0 substantially overproduce the number of clusters, so the 
consideration of these models here is essentially motivated by didactic purposes, as an example 
of pure power law spectra. Affecting a different normalization - for example, normalizing 
directly to the mass or temperature distributions - can, in principle, yield different results, due 
to the different implied shapes and amplitudes of N(M ) z). 

A comparison of the number counts for flat and low density CDM cosmologies shows that the 

counts increase with decreasing VL$. With our normalization to the XRLF, lowering f2 has the 

effect of lowering N(M, z) (which is proportional p) and reducing the cutoff mass. As a result, 

the local abundance of clusters in a low density universe is lower than in the Q = 1 case, if 

we assume the same values for h, X and b5 v . This effect is counterbalanced by lowering bS v and 

increasing h and/or X (see Table 5). In addition, in a low density universe N(M, z) evolves less 

than in a flat universe (see Fig.5). We stress that our scaling laws imply that the core radius 
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is larger in low density a universe as compared to a flat model. With all the other parameters 
fixed, the net flux from the cluster is then increased. Altogether, without gas evolution, we 
do expect a factor m 5 more clusters in a low density than in a critical CDM universe. This 
difference is reduced when including evolution of the IC gas which has the effect of making the 
population of predicted clusters a more local (z ^ 0.4) one. Note the strong dependence of the 
counts on the assumed value of the Hubble constant. 

Among the various assumptions and simplifications made in our treatment, the uniform, 
spherical cluster collapse - as described in the Press & Schechter (1974) formalism - is a major 
one. Another major uncertainty is the nature and degree of IC gas evolution. The two limiting 
cases considered here - no evolution, and maximum degree of evolution which is still consistent 
with the EMSS distant XRLFs - are likely to span a reasonably realistic range. It should 
also be noted that our scaling of the gas fraction to a value of 0.1 in a rich, local cluster may 
be conservatively low, judging by some observational indications of a value of up to 30% (see 
White et al. 1993). These uncertainties are not expected to affect our conclusion that the sky- 
averaged optical depth for electron scattering in the hot IC gas is quite low in all the models 
considered here. Obviously, the range of predicted values of r constitutes the minimum level 
of scattering of the CMB eversince the epoch of last scattering. Reionization before cluster 
formation may have resulted in a much higher level of scattering (Tegmark & Silk 1995). 

Another simplification made in our treatment is the use of the non-relativistic expression for 
the S-Z effect. As we have noted already, in the exact relativistic treatment the dependence of 
the S-Z effect is not linear in y, and the the deviation of the intensity change from the non- 
relativistic value is appreciable at high frequencies and gas temperatures (Rephaeli 1995b). The 
major contribution to the integrated effect due to the full cluster population comes from low 
temperature clusters. The full relativistic calculation can be affected by means of a frequency- 
dependent correction factor (applied to the non-relativistic expression for the intensity change) 
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which can be calculated at some mean, population-weighted value of the gas temperature. 
Since we expect this mean to be close to the lower end of the observed range, 3 15 keV, the 
correction will generally be small except at very high frequencies. For example, if this mean 
temperature is ~ 5 keV, then the correction factor is a few percent in the Rayleigh- Jeans part 
of the spectrum, and ^ 20% at the 200 240 and 600 GHz ranges. The exact correction 
factors can be calculated at each frequency and gas temperature using the expressions given in 
Rephaeli (1995b). 

The cluster number counts in the mm and sub-mm regions are of interest for upcoming high- 
sensitivity CMB anisotropy experimentsCMB anisotropy experiments such as BOOMERANG 
and C/S. The number counts shown in Table 3 are optimistic in several ways. First, sky confu- 
sion was not taken into account; its inclusion will obviously reduce our predicted numbers in a 
way which will largely depend on the degree of sensitivity in modeling emission from Galactic 
dust and far-IR emission from other galaxies. The Galactic disk region will reduce useful sky 
coverage to ~ 80%. The integrated emission from galaxies is not known, but in some models 
with strong luminosity or density evolution the predicted intensity levels (e.g., Beichman and 
Helou 1991, Wright et al. 1994) may well exceed those corresponding to Comptonization . A 
quantitative comparison is not warranted at this stage because of the high degree of uncertainty 
in these models. In the analysis of actual data it will likely be possible to separate out the S-Z 
component based on its unique spectral shape and its larger characteristic spatial scales. For 
now, we include in Table 4 the numbers of clusters with AF U > 3F^, mse . From these we can 
predict that if sky-confusion is minimal, then a full sky map - which can be generated by C/S 
at 400 GHz (i.e. , with a resolution of 4'. 5) after one year of operation - should include as many 
as ~ 10 3 clusters detected at the 3a statistical significance level. 
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Appendix 



The expansion rate of a low density universe is given by the Friedmann equations: 

{ H 2 (l - n )R ; A = ' 

> 

[ H 2 {l - {lo)R 2 ; A = l - fi 



R2= 8nGp R2 + 



(25) 



R 2 = —PR 2 + 



(26) 



The second terms in the rhs of the previous expressions describe the curvature and the cosmo- 
logical constant contributions. 

Let us consider a spherical, homogeneous perturbation. If bound, this perturbation evolves 
according to the following relations: 

-c 2 ; A = 

1 -Q )H 2 R 2 -c 2 ; A=l-Q 

Here p is the perturbation overdensity and the curvature terms (oc c 2 ) are positive. Note 
the repulsive effect of a non vanishing cosmological constant (oc 1 — fi ) on the perturbation 
evolution. 

Let us assume that at some initial time U the fluctuation has the same size and expansion 
rate of the background. By subtracting Equation (27) from Equation (28) we find the conditions 

f -H 2 (l - n )R 2 ; A = 0) 

c 2 = ^ P AR 2 + (27) 
6 { 0; A = l-ll 0; 

Here <5j and Ri are the density fluctuation and proper size of the perturbation, is the back- 
ground density at U and R m is the proper size of the perturbation at turnaround. 



Substituting Equation (29) in Equation (28) we obtain 



3(1 — Q )HqRRI; 



3RR 2 = 8irGp~R - ^GpARfR + { 



A = 



k 3(1 - n )H 2 R 3 ; A = 1 - Q 



(28) 
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Imposing the turnaround condition, R = 0, one finds in both cases 

A = 



f2o r Rq 



R m - < 



l — fio Ri 



V 

{ 2 + 7] 



; A = 1 - fi 



(29) 



where r\ = 3(1 — Q )HQ/(4irGp m ). Substituting Equation (31) in Equation (30) yields: 



r 2 



3RR 2 = 3H 2 (1 - tt )(R m -R)x 



• d2 . 
Jn m , 



A = 



—R 2 — R m R H — Rm)'i A — 1 — fio 



(30) 



When integrated, the previous expressions read: 

x' 



H Jl-tt t= I 



rj/2 £ dx' 

rx 

/ dx' 
Jo 



(1 - X') 



A = 



x 



1/2 



; A = 1 - fl 



(31) 



.(1 -x')(2/r]-x' -x 12 ). 
where x = R(t)/R m . Analytical integration of Equation (33) with A = provides the standard 
cycloid solution. 

At U, by definition Xi = R(U)/R(t m ) <C 1: this allows a Taylor expansion of the integrands 
of the rhs of the previous equations and an analytical estimate of the integrals. One gets: 



Rrn\ 3/2 = 2 frf 1 

Ri) ~ 3V2# 0V /i -n u 



1 + > 



10 V3 V 2 H oy /l - n ti 



-2/3-1 



x < 



10 V 2A3V2 J fWl-fio*i 



-2/3- 



A = 

A = 1 - fi 

(32) 



We can now write the initial overdensity of the fluctuation: 

3/2 frf 1 N- 2 / 3 



R ^ 3 
Ri 



1 + 



5V3V 2H Q ^T^TT U 



6nGt 



x < 



5V 2/131/2^1-^, 
18 



-2/3 



A = 



; A = 1 - n 



(33) 



As both in the open and flat case pi = l/(6nGt?), given that pi = p;(l + Si), one finally gets: 

f 1; A = 1 



3/3 2 



5 V2 V v 



Hoy 1 — &ot, 



2/3 



X < 



1 + 



7/ 



A = 1 - fin 



(34) 



The previous expressions give the initial density contrast (at some arbitrary early time tj) 
of a perturbation that will turnaround at a given t m . This information is contained in the 
r) parameter through Equation (33). Substituting t — t m in eq. (31) which implies x — 1, 
we obtain a relation between the turnaround time and r\. Given the parameter r\ we have 
(from its definition) the value p m of the perturbation overdensity at the turnaround time. We 
assume that the perturbation eventually virializes at t v = 3t m and in this regime reaches a final 
dimension R v = R m /2 in the open case, and R v = R m (l — r//2)/(2 — rj/2) in the A = 1 — f2 
case (see Lahav et al. 1991). Having determined the virialization time, t v , and the background 
density at that time, we can finally evaluate the non-linear overdensity of a virialized structure 



P(tv) 

In an open universe we can write the following analytic expression: 

18tt 2 1 



(35) 



" n Hffi (i + Zv f ' 

where t v and related through standard time-redshift relations. 



(36) 
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Figure captions 



Figure 1. The dependence of y from the maximum redshift, z max , of integration [cf. Equation 
(16)]. A minimum mass cutoff M min = 10 125 h _1 M has been chosen, corresponding to the 
smallest systems we consider here. Continuous and dotted curves refer to cases with and without 
IC gas evolution. 

Figure 2. The integration domain of iV(> AF V ) is shown for the case of a standard CDM 
model (fi = I, h = 0.5, n = 1) and for two different channels of the C/S experiment. The 
dotted contours in the LogM-Logz plane represent isocontours of the quantity N(M, z)-M-z and 
are spaced by a factor of 10, the outer contour corresponding to unity. The region contributing 
to the source counts is to the right of the lines AF U = 161 and 19.7 mJy for the C/S's 140 and 
400 Ghz channels, respectively. Heavy continuous and dashed lines refer to the cases with and 
without IC gas evolution. 

Figure 3. Cluster number counts for the C/S 400 GHz channel for different angular resolu- 
tions (4.5 arcmin: continuos lines; 10 arcmin: dotted lines; 30 arcmin: short-dashed lines; 60 
arcmin: long-dashed lines) in different models with (Figure 3a) and without (Figure 3b) IC gas 
evolution. 

Figure 4. Same as if Figure 3 but for the C/S 140 GHz channel. In this figure we do not 
plot the predictions for the 4.5 arcmin resolution. 

Figure 5. Comparison of the cluster number counts predicted in different models of structure 
formation (MDM: continuous line; open CDM: dotted line; low density, vacuum dominated 
CDM: short-dashed lines; standard CDM: long-dashed lines) for the C/S 400 GHz channel 
with (panel a) and without (panel b) IC gas evolution. 

22 



Figure 6. The redshift distribution N(> z) of the cluster number counts in different models 
(MDM: continuous line; open CDM: dotted line; low density, vacuum dominated CDM: short- 
dashed lines; standard CDM: long-dashed lines) for different sensitivities: AF U = fO and 20 
mJy (Figure 6a) and AF U = 30 and 60 mJy (Figure 6b) for the C/S 400 GHz (a FWHM = 4.5 
arcmin) and 140 GHz ((Jfwhm = 10 arcmin) channels. In each figure panel a) and b) refer to 
predictions without IC gas evolution. 

Figure 7. The distribution jV(> 9fwhm) of the typical angular dimensions of the clusters 
above the limiting flux for the C/S 400 (Figure 7a) and 140 (Figure 7b) GHz channels. Thin 
lines, with (continuous) and without (dotted) IC gas evolution, refer to the distribution of 
intrinsic dimension. Thick lines, with (continuous) and without (dotted) IC gas evolution, refer 
to 4.5 and 10 arcmin smoothing of the C/S maps. 
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Table 2 



Models 


<y> in units of 10" 6 




with ICM evolution 


without ICM evolution 




10 13 M o 


10 13 - 5 M 


1 A 

10 14 M o 


10 13 M o 


IOC 

10 13 ' 5 M 


1 A 

10 14 M o 


CDM 














£2=1.0; h=0.5 














n=l 


0.34 


0.30 


0.23 


1.17 


0.86 


0.52 


n=0.8 


0.30 


0.27 


r\ fit 

0.21 


0.92 


0.71 


0.45 


CDM+A 














£2=0.2; n=l 
h=0.5 
h=0.6 
h=0.7 


0.27 
0.22 
0.18 


0.26 
0.20 
0.16 


0.23 
0.17 
0.13 


0.99 
0.88 
0.77 


0.82 
0.70 
0.58 


0.61 
0.48 
0.37 


CDM open 














Q=0.2;n=l 














h=0.5 


0.20 


0.18 


0.16 


1.16 


0.89 


0.60 


h=0.6 


0.14 


0.13 


0.11 


0.88 


0.63 


0.38 


h=0.7 


0.12 


0.11 


0.08 


0.77 


0.52 


0.28 


MDM 














Q v =0.3;h=0.5 














n=1.0 


0.07 


0.06 


0.05 


0.15 


0.12 


0.08 


n=1.2 


0.10 


0.09 


0.07 


0.24 


0.19 


0.12 


n=1.4 


0.16 


0.14 


0.11 


0.41 


0.31 


0.20 


BDM 














£2=0.1; n =-l 














h=0.5 


0.43 


0.41 


0.37 


27.3 


15.0 


7.25 


h=0.6 


0.33 


0.31 


0.27 


18.9 


9.76 


4.27 


h=0.7 


0.25 


0.23 


0.19 


13.2 


6.41 


2.54 


BDM 














a=0.4; n=-l 














h=0.5 


0.46 


0.42 


0.34 


5.19 


3.04 


1.53 


h=0.6 


0.35 


0.31 


0.23 


3.79 


2.08 


0.94 


h=0.7 


0.32 


0.27 


0.19 


3.38 


1.78 


0.74 



Table 3 



Models 


<x> in units of 10" 4 




with ICM evolution 


without ICM evolution 




10 13 M o 


10 13 - 5 M 


10 14 M o 


10 13 M o 


10 13 ' 5 M 


10 14 M o 


CDM 














£2=1.0; h=0.5 














n=l 


1.53 


1.04 


0.56 


6.58 


3.34 


1.34 


n=0.8 


1.27 


0.88 


0.49 


5.05 


2.68 


1.14 


CDM+A 














Q=0.2; n=l 
h=0.5 
h=0.6 
h=0.7 


0.84 
0.79 
0.75 


0.65 
0.59 
0.53 


0.45 
0.37 
0.30 


4.75 
4.56 
4.28 


2.74 
2.47 
2.18 


1.38 
1.13 
0.90 


CDM open 














Q=0.2;n=l 














h=0.5 


0.57 


0.43 


0.29 


4.72 


2.50 


1.14 


h=0.6 


0.50 


0.36 


0.21 


4.01 


1.96 


0.79 


h=0.7 


0.48 


0.32 


0.17 


3.70 


1.69 


0.61 


MDM 














Q v =0.3;h=0.5 














n=1.0 


0.34 


0.22 


0.12 


0.91 


0.50 


0.22 


n=1.2 


0.50 


0.33 


0.17 


1.49 


0.78 


0.33 


n=1.4 


0.76 


0.50 


0.26 


2.51 


1.29 


0.53 


BDM 














£2=0.1; n=-l 














h=0.5 


0.81 


0.69 


0.51 


58.3 


23.8 


8.42 


h=0.6 


0.77 


0.62 


0.43 


45.1 


17.2 


5.55 


h=0.7 


0.69 


0.53 


0.33 


35.0 


12.5 


3.63 


BDM 














£2=0.4; n=-l 














h=0.5 


1.63 


1.21 


0.73 


22.2 


9.31 


3.28 


h=0.6 


1.45 


1.01 


0.54 


17.8 


6.93 


2.15 


h=0.7 


1.38 


0.92 


0.45 


16.2 


5.98 


1.69 



Table 4 



Models 


N (>F v noise ) 




with ICM evolution 


without ICM evolution 




4.5' 


5' 


10.5' 


4.5' 


5' 


10.5' 




400GHz 


150GHz 


140GHz 


400GHz 


150GHz 


140GHz 


CDM 














£1=1.0; h=0.5 














n=l 


4981 


7 


146 


18258 


25 


286 


n=0.8 


DZUO 


o 
o 


1 81 
lol 




ZD 




CDM+A 














£1=0.2; n=l 
h=0.5 
h=0.6 
h=0.7 


7U94 
4698 
3194 


11 
7 
5 


350 
200 
118 


26336 
19722 
14583 


38 
28 
20 


611 

376 
233 


CDM open 














Q=0.2;n=l 














h=0.5 


5518 


9 


260 


32821 


47 


631 


h=0.6 


3065 




121 


19426 


27 


297 


h=0.7 


2157 


3 


75 


14735 


20 


195 


MDM 














Q v =0.3;h=0.5 














n=1.0 


1348 


2 


65 


3176 


5 


101 


n=1.2 


1745 


3 


74 


4526 


7 


122 


n=1.4 


2599 


4 


101 


7422 


11 


175 


BDM 














£2=0.1; n =-l 














h=0.5 


12114 


19 


494 


414779 


511 


2652 


h=0.6 


7511 


7 


273 


237724 


284 


1295 


h=0.7 


4525 


7 


150 


138250 


161 


666 


BDM 














Q=0.4; n=-l 














h=0.5 


6930 


10 


226 


26597 


60 


525 


h=0.6 


4145 


6 


120 


29466 


37 


302 


h=0.7 


3611 


5 


104 


27995 


35 


283 



Table 5 



Models 


N (> 2F v noise ) at 400 GHz 




with ICM 
evolution 


without ICM 
evolution 


CDM 

£2=1.0; h=0.5 
n=l 
n=0.8 


1063 
1220 


3459 
3792 


CDM+A 

£2=0.2; n=l 
h=0.5 
h=0.6 
h=0.7 


2068 
1259 
786 


6828 1 
4648 

3111 j 


CDM open 

£2=0.2;n=l 

h=0.5 
h=0.6 
h=0.7 


1629 

810 

528 


8872 
4656 
3261 


MDM 

£2 v =0.3;h=0.5 
n=1.0 
n=1.2 
n=1.4 


364 
441 
633 


771 

1027 ■ 
1620 


BDM 

£2=0.1; n=-l 

h=0.5 
h=0.6 
h=0.7 


3449 
1975 
1105 


95100 

49832 1 
26549 


BDM 

a=0.4; n=-l 

h=0.5 
h=0.6 
h=0.7 


1600 

875 

750 


9011 

5205 ! 
4889 



Table 1 



1V1UU-C1S 


Fit parameters 




h 


I 


bd v 


n 




TDM 0-1 n 
n=l 


0.5 


1 


2.5 


0.2 


1.2 


n=0.8 


0.5 


1 


2.45 


0.2 


1.2 


LUM+A; £2=0.2; n=l 


0.5 


1 


1.3 


0.2 


1.45 




0.6 


1 


1.4 


0.2 


1.6 




0.7 


1 


1.5 


0.2 


1.7 


fDM- O-D ?• n— 1 


0.5 


1 


1.4 


0.2 


1.6 




0.6 


1 


1.55 


0.2 


1.7 




U./ 


i 


l.OO 


U.z 


1 Q 

l.O I 


MDM; Q v =0.3 

n=1.0 


0.5 


1 


2.8 


0.2 


0.8 


n=1.2 


0.5 


1 


2.75 


0.2 


0.9 


n=1.4 


0.5 


1 


2.65 


0.2 


1.0 


BDM; a=0.1; n=-l 


0.5 


1 


0.7 


0.2 


2.0 




0.6 


1 


0.8 


0.2 


2.1 




0.7 


1 


0.9 


0.2 


2.25 


BDM; a=0.4; n=-l 


0.5 


1 


1.5 


0.2 


1.7 




0.6 


1 


1.7 


0.2 


1.85 




0.7 


1 


1.8 


0.2 


2.0 



